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Abstract 

We propose a method for construction of Darboux transformations, which is a new devel¬ 
opment of the dressing method for Lax operators invariant under a reduction group. We apply 
the method to the vector sine-Gordon equation and derive its Backlund transformations. We 
show that there is a new Lax operator canonically associated with our Darboux transformation 
resulting an evolutionary differential-difference system on a sphere. The latter is a generalised 
symmetry for the chain of Backlund transformations. Using the re-factorisation approach and 
the Bianchi permutability of the Darboux transformations we derive new vector Yang-Baxter 
map and integrable discrete vector sine-Gordon equation on a sphere. 


1 Introduction 

Lax integrable partial differential equations have natural connections with integrable differential- 
difference and discrete equations via Darboux transformations. Namely, given a Lax representation 
for a partial differential equation, we can systematically construct Darboux transformations whose 
Bianchi permutability condition leads to an integrable difference equation, while the corresponding 
Backlund transformations are (often nonlocal) symmetries of these difference equations and are 
integrable differential-difference equations in their own right. 

In this paper we make steps towards the development of a systematic approach to construction 
of Darboux transformations, associated difference systems and integrable maps based on a natural 
extension of the dressing method. The motivation for this line of research is to cooperate the 
reduction groups mmm of Lax representations of integrable partial differential equations to 
integrable difference equations and corresponding Darboux transformations. The aim of our project 
is to describe all elementary Darboux transformations for Lax operators on Kac-Moody algebras 
and automorphic Lie algebras with finite reduction groups. Recently, the authors of HE] have 
completed a comprehensive study for the Lax operators of the nonlinear Schrodinger equation type 
and they derived some new discrete equations and new Yang-Baxter maps. In the paper [6| the 
method was applied to the case of Lax operators with the reduction group generated by inner and 
outer automorphisms of the Lax representation for two dimensional Volterra chain. 

In this paper we study differential-difference and partial difference equations associated with 
the vector sine-Gordon equation 
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where dependent variables o: = (a^, • • • , a^)'^ is n-dimensional real vector and /3 G M. Here and in 
what follows the upper index T denotes the transposition of a vector or a matrix and we denote 
= Ictp. Sometimes it is convenient to write /3 and a as an n + 1-dimensional unit vector 
v'^ = ( /I, ). Thus the vector sine-Gordon equation is an integrable partial differential system 

of equations with the vector dependent variable is on a sphere v £ S"' = {x £ | x'^x = 1}. 

Equation © has a long history, it first appeared in [7] viewed as a reduction of the two- 
dimensional (T-model [8j. Its Lax representation was given in [9] and its Lagrangian formulation in 
m- Later, this equation reappeared in the study of moving frames in the Riemannian geometry 
m- The dressing method for construction of multi-soliton solutions of © and study of soliton 
interactions has been recently developed in m- 

In this paper we use the Lax representation (fTUD , (fTTD proposed in m and studied in [T2] . For 
these Lax operators we derive a Darboux transformation (Theorem [1]) with Darboux matrix 

'illy 'Liy 1 I 

M„{u) = In +2 +-. - —Qa){a-——a){aQ, {a = {i,u^), u = - -^—r, (2) 

A — IV \ + IV 1171 + 

where In +2 is the identity matrix of size n + 2, Q = diag(—1,1,..., I), and rigorously prove that 


Dt 


f Vi+V 

Vl^i +1’ 


/3i 13 




2v 


+ V 


(qi -I- q) , 


(3) 

(4) 


is the Backlund transformation which depends on an arbitrary constant parameter 0 yi: G M and 
relates two solutions v and vi of the vector sine-Gordon equation ([1]). Equations (l3|) and © can 
also be seen as two non-evolutionary differentiaEdifference equations on a lattice with variables 
Vj £ 5”, j £11 and with a shift operator Sp : Vj i—)• 'Uj+i, assuming vq = v. 

For the Darboux matrix ([2]) there exists a new Lax operator 


_ _ . Qa-i){a - a){a-iQ . a-i){aQ - Qa){a-i 

{X — iv){aQa^i) {X + iv){aQa-i) ’ 


«-i) =‘5j,Ha)) 


(5) 


canonically associated with it. The existence of a Lax operator which is canonically associated 
with a Darboux matrix is quite remarkable. We hrst observed it in [6] where a local generalised 
symmetry for non-evolutionary differential-difference equations (Backlund transformations for two- 
dimensional Volterra chain) has been found. Our analysis of the variety of Darboux-Lax integrable 
differential-difference systems shows that local symmetries often correspond to such operators. 
They deserve more attention and we are planning to develop this direction of research in future 
publications. 

The compatibility of the Darboux matrix ([2]) and Lax operator ([5]) also results in an evolutionary 
differential-difference equation 


Dr{v) 


I’U-l + v\'^{vi -|- — |l>l -|- -|- v) 

{vi + v'^){v_i ^») -b |^’-l + 


Vk = s^{v), 


( 6 ) 


which is a local symmetry of the nonevolutionary integrable systems ([3]) and (jH). System ([3D is 
known, and it was found by Adler m in his classification of isotropic integrable Volterra-type 
lattices on the sphere with generalised symmetries. In this paper we equip this system with a Lax- 
Darboux representation and connect it to the vector sine-Gordon equation ([ID and its Backlund 
chains ([3D and ([ID. 
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The Bianchi permutability condition for two Darboux transformations with distinct parameters 
fj. 7 ^ ± 1 / resulting in two shift operators S^,, leads to the integrable discrete equation 


' 1 ^ 1,0 = -V + 2X(X'^, v) 


(7) 


where = Vn,m and 

^ + 2fi{u + /r(x^y))y 

1/2 + _g 2 /iz/(x^y) 


— - - - — ? y “ — - - 

1^1,1+^’ 0 , 1 !’ |^^0,l+^’ 


Equation ([7|) can also be uniquely resolved with respect to the variable uoq. Backlund transforma¬ 
tions m, 0 and evolutionary equation ([6]) are its non-local and local symmetries. 

Moreover, the re-factorisation of the product of two Darboux matrices M^(x.)M^(y) = M^{Y)Mi,(X) 
leads to a new involutive Yang-Baxter map Ru,^i ■ x x 


R 





X 

Y 


/ (t/2-^2)x + 2/r(z/ + Mx^y))y \ 
1/2 -h ^2 _g 2 /iz/(x'^y) 

(^2-^2)y + 2^(^ + ^(xTy))x 
\ 1/2 -h /u2 + 2^z/(x'^y) / 


x,y,X,Ye5T (8) 


The arrangement of this paper is as follows: In Section [21 we give basic definitions related 
to the Lax- Darboux scheme such as Darboux transformation, the shift operators and Backlund 
transformation. In Section [3l we derive the Backlund transformations for the vector sine-Gordon 
equation ([T|) together with the Darboux matrix. In Section jH we construct a local generalised 
symmetry for these non-evolutionary equations using the Darboux matrix. Meanwhile, we provide 
the Lax representation for this symmetry flow, which is one of isotropic integrable Volterra-type 
lattices on the sphere m- In Section El we use the re-factorisation of two Darboux matrices to 
derive the corresponding Yang-Baxter map and integrable discrete vector sine-Gordon equation, 
which are new to the best of our knowledge. 


2 The Lax-Darboux scheme 


In this section, we recall the Lax representation of the vector sine-Gordon equation and introduce 
some basic definitions such as Darboux transformation and Backlund transformation. 

The vector sine-Gordon equation ([1]) is equivalent to the compatibility condition [T, M] = 0 for 
two linear problems m, m 

= 0, MT = 0, (9) 

where 


C = D^- XJ -U and A = Dt + A"V, 


/ 0 1 

J= -I 0 0^ 

V 0 0 On 


/ 0 0 

[/ = 0 0 

\ 0 a^//3 On 


V = 


0 v'^ \ 

-V On+l ) ’ 


( 10 ) 


( 11 ) 
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where 0 is n-dimensional zero column vector and 0^ denotes the kxk zero matrix. Without causing 


confusion, we sometimes simply write 0 instead for either of them. 

A Darboux transformation is a linear map acting on a fundamental solution 

= M'h, det M / 0 (12) 

such that the matrix function 'h is a fundamental solution of the linear problems 

= 0, A'^ = 0 (13) 

with new “potentials” a and (5. The matrix M is often called the Darboux matrix. From the 
compatibility of (fT^ and (fT^ it follows that 

DAM) = (AJ + U)M - M{XJ + U) = X[J, M] + UM - MU] (14) 

Dt{M) = {-X-^V)M - M(-A"V) = -A“^(FM - MV), (15) 


Equations lead to Backlund transformations for the vector sine-Gordon equation, which relate two 
solutions a, f3 and a, (3 of ©• 

A Darboux transformation maps one compatible system ([9]) into another one m- It defines a 
map S : V >-^v. The map ()12p is invertible (det M 7 ^ 0) and it can be iterated 

T S T S <5 S 

••• . 


We introduce notations 


... T-i = T, To = T, Ti = T, T2 = T,..., 

. . . U_1 = V, Vq = V, Vi =V, V 2 =v, . . . . 

In these notations the maps S and 5“^ increase and decrease the subscript index by one, and 
therefore we shall call it a <S-shift, or shift operator <S. The resulting Backlund transformations 
from the Lax-Darboux representations (1141) and (1151) are integrable differential difference equations. 

A Darboux transformation with a parameter ^ denoted by M^ results in the shift. If 
we also consider a Darboux transformation with a different choice of the parameter then the 
corresponding shift we denote S^,. Commuting shifts act on lattice where with the vertex (n, m) 
we associate the variable Vn^m = Commutativity of the shifts is (Bianchi permutability) 

SAM^M^ - = 0 (16) 

is equivalent to a system of partial-difference equations and we call it Darboux representation 
for this system of partial-difference equations. In literature such representation and Darboux 
matrices sometimes referred as discrete Lax representation and discrete Lax operators respectively. 
Differential difference equations (|14ll and (|15p are the symmetries of this partial-difference equation. 

We can also consider the re-factorisation of a product of two Darboux matrices, which leads 
to the Yang-Baxter map. In next sections, we construct the Darboux matrices for the vector 
sine-Gordon equation ([T]), and further produce its Backlund transformations, the associated Yang- 
Baxter map and the integrable difference equation. 
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3 Invariant Darboux matrix under the group Z2 x Z2 x Z2 


The Lax operators L and A are invariant under the reduction group Z 2 x Z 2 x Z 2 generated by 
three transformations i r and s satisfying = id. In particular, for operator L these 

transformations are: 


1 : L{\) -0{\) = /:(A), 

r : £(A) ^ £*(A*) = £(A), 
s:£(A)^Q£(-A)Q = £(A), 


(17) 

(18) 
(19) 


where T^(A) is the formally adjoint operator defined by T^(A) = — AJ^ — matrix Q = 

diag(—1,1,... , 1) and * denotes the complex conjugation. 

We assume that Darboux matrix M(A) is a rational function of the spectral parameter A and 
is invariant with respect to the reduction group CZD-CID action, namely. 


1 : M(A) ^ = M(A), 

r : M(A) ^ = M(A), 

s : M(A) ^ QM{-X)Q = M(A). 


( 20 ) 

( 21 ) 

( 22 ) 


Notice that the action of the automorphism 1 m is different from (HZl) since M is an element of 
a Lie group rather than Lie algebra. 

It is easy to show (see Proposition 1 in [12]) that a A-independent Darboux matrix results in a 
constant linear map. 


Proposition 1. [T2| Assume Darboux matrix M is invariant under the group Z 2 x Z 2 x Z 2 and it 
is independent of spectral parameter A. Then M is a constant matrix and of the form 


M = ± 



(23) 


where 0 is n-dimensional zero column vector and a constant (x, t - independent) matrix D € 0(re, M). 
This implies that the constant Darboux matrix defines the transformation 


jd = fi and a = Da. 


(24) 


corresponding to 0(n,M) point symmetry of the vector sine-Gordon equation ([T]). 

A rational in A and reduction group invariant (I20p - (I22D matrices M must have poles on the 
orbits of the reduction group. They can be represented as a product of matrices which have poles in 
a single orbit. Here we shall study reduction group invariant Darboux matrix M,^ with two simple 
poles at A = dziiz, n ^ 0, € R. It is the minimal possible number of simple poles for a reduction 

group invariant Darboux matrix and it corresponds to the dressing matrix for a single kink solution 
j. It is of the following form (cf. Proposition 2 in |12jL 


Mu — In+2 + 


IV 


\ — iu 


-A- 


liz 


X + in 


—QAQ, A A 0. 


(25) 
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Proposition 2. If matrix is invariant under reduction group generated by i, r and s, it can 
be represent as 


M^{u) =In+ 2 + , . Qa){a- . a){aQ, {a = {i,u^), (26) 

X — IV X + iv 

where u is a unit length real vector u £ S^. 

Proof. Note that the given matrix is invariant under the action of s. It immediately follows from 


(120)) that MyM^f = In+ 2 , which implies that 

AA^ = A^A = 0; (27) 

{In+2 ~ + ^(-In+2 “ Q) = 0. (28) 

These identities correspond to the vanishing of the second and first order poles at A = hiv, 
respectively. The invariance under the action of r implies 

A* = QAQ, (29) 


So the matrix A satisfies the conditions of Lemma [T] in Appendix. Hence, it is of the form 


A = h( * 

Y m 


i u ) , 


0 7 ^ /i € M, m, u € |m| = |u| = 1. 


Substituting it into (l28)) . we obtain 

/i = 1, m = u. 

Thus we have A = Qa){a, from which the result in the statement immediately follows. □ 


This proposition confirms that the kink solutions of rank 1 for the vector sine-Gordon equation 
© obtained in [12] are indeed generic and there are no kink solutions of higher ranks. 

Proposition 3. Matrix M^{u) \2b\) satisfies the compatibility condition (ESP if and only if 



VM^,{u) = M^{u)V, 

(30) 


1 _ 1 

(31) 


Dt {Qa){a) + — VQa){a - —Qa){aV = 0, 

IV IV 

where 

M'^{u) = lim M^{u). 

A^O 



Proof. Taking the residue at A = 0 on both sides of (fT^ . we obtain (1301) . Taking the residue at 
X = iv on both sides of (fTCD . we obtain (I^T]1 . The residue at A = —iv will also vanishes due to the 
reduction group. □ 


We now convert the above conditions to the conditions for the components of vector (a = {i, u^). 
Proposition 4. Matrix M^{u) given by satisfies [W\) if 


u = 


± 


V + V 

V + V 


(32) 
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(33) 


Proof. It follows from (I26p that 

M°(it) = In +2 - Qa){a - a){aQ = (M°(w))'^ = (M°(u))“^ 
Thus we have V = M^{u)VM^{u). Using (fTTll . we get 

V + V = 2{u'^v)u. 

It follows that + 1) = ||u + up. Thus we have u = ±|=i^. 


□ 


In the rest of the paper we assume that v -\-v ^ 0 and choose the positive sign in (13211 . We now 
compute the Backlund transformation and summarise the result for the Darboux matrix. 


Theorem 1. The vector sine-Gordon equation ^f) possesses a Backlund transformation 




=-(u — u); V = 


\v V 

where v is constant. The corresponding Darboux matrix is 


a 

OL 


^ V , IV ^ , '• T\ V -\-V 


M^{u) = In +2 +-. - —Qa){a— . a){aQ, {a = {i,u ), u = — 


X — iv X + iv ' ' ' V - |.j; _j_ .y| 

Proof. We derive the Backlund transformation from ()31ll . which is equivalent to 


1 


lU 


Q { at) —Q VQa) {a + Qa) {at - —{aV = 0. 


1 


Therefore, there exists a scalar 7 such that 

at) + ^Q~^VQa) = 0)7; 

IV 

Taking the sum of these two equations we get 

I 


IV 


at) + —V a) = -0)7. 

IV 


at) + — (Q-VQ + U) a) = 0. 


(34) 


(35) 


(36) 


(37) 


that is, 


(u^ — v^)u = 0; 


Dtu =-(u — u). 

2v^ ’ 

The hrst equation is satisfied due to ([52]) and |u| = |ui| = 1. The second one is the required 
equation ([Mil . It is easy to check that equations (l36[) and (l3711 (with 7 = ^v'^u) are satisfied due 
to (IMll. □ 


Using the compatibility condition ([14p between the Darboux matrix (13511 and matrix U, we get 
the following Backlund transformation with respect to x: 

Corollary 1. For the Darboux matrix Mjj{u) given in Theorem{J\ the compatibility condition 
leads to the following Backlund transformation (with respect to x): 


Dxja) 


Dxja) 

P 


2v 


\v + v\ 


( Q + q) 


(38) 
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Proof. Substituting the Darboux matrix My{u) into (1141) and comparing the residues of all poles 
from both sides, we have 


U = U + iv [Qa) {a — a) {aQ, J] ; 

Dx (Qa) {a) = UQa) {a — Qa) {all + iu [j, Qa) (a] 


By direct calculation, we obtain the identity (138)1 in the statement from (l39l) . 

We now check that ini) satisfies under the condition of (|38p . We rewrite it as 


Q [Dx{a)) — Ua) + iuJa)) {a + Qa) {{Dx{a) + {all + iu{aJ) = 0. 
Differentiating + Dp with respect to x and using I3x = —^cc^ax, we have 


Dx(\v + v\) = 


2u 


1 


(d'^ + )Dx{v + v) = 


T\ n 


\v + V 

t-T -5„T 


1 , p—T 15 T\ ( OLx 

- dpOL —PCX ) —-=■ 

l^ + ^l (3 


{ficx^ — Pcx ^) (q + a) = v{j3 — (3). 


Using it and 


|d + Dp 

, by direct calculation we obtain 


(39) 

(40) 


(41) 


{Dx{a) + {all + ii'{aJ = — 7 =-—^ {a; 

\v + d| 

Dx{o)) — Ua) + ivJa) = ^ a). 

|d + d| 

Substituting it into (I4T|) . we see the identity is valid and thus we proved the statement. □ 

When n = 1, we take a = sin6( and f3 = cos 9. The vector sine-Gordon equation ([T]) reduces to 
the scalar sine-Gordon equation 


Oxt = sin 9 


(42) 


We show that (f34|l and ([381) give 
In this case, we have j=^ = 
identities, we obtain 


us its the well-known Backlund transformation. 

e+e \ 

. Substituting it into ([34)1 and using trigonometric 



9t + 9t = - sin 


e-e 

2 • 


Substituting the above a and (3 into the Backlund transformation (13811 . we get 


(43) 


9x — 9x = -- 


2 d 


2 cos 


^—(sin0 -|- sin0) = —2d sin 


e+e 


e-e 


(44) 


Backlund transformations (13811 . (14411 provide us with integrable discretisations of the vector 
m and scalar (j42p sine-Gordon equations. There exists another integrable discretisation of these 
equations, recently published in [T4], which has a different nature and is not directly related to the 
Backlund transformations of the vector or scalar sine-Gordon equations. 









4 A local symmetry of the Backlund transformation 


As discussed in Section [2] the Backlund transformations and (1551) can be viewed as integrable 
differential-difference equations. In this section, we show how to construct one of their generalised 
symmetries using the Darboux matrix (|35p . This symmetry is an integrable evolutionary equation, 
which belongs to the list of integrable Volterra-type lattice systems on a sphere satisfying the inte- 
grability conditions following from the existence of a formal recursion operator m- We construct 
a Lax operator, canonically associated with the Darboux matrix (|35l) which forms Lax-Darboux 
representation for this local symmetry. 

The Lax operators A, L (fTOl) have simple poles at points of the degenerate orbits A = 0, oo 
and are invariant with respect to the reduction group generated by the transformations (|17l) - (ll9p . 
Starting from the Darboux matrix M^{u) (1351) we construct a new Lax operator 

B^ = Dr-U^, (45) 


which is invariant with respect to the reduction group and matrix Uy has the same set of poles as 
the logarithmic Frechet derivative of the Darboux matrix. Thus, in the case of the Darboux matrix 
My{u) (|35p matrix lAi, has the form 


U,{X) 


ivB iuQBQ 
\ — iv \ + iu 


(46) 


and satishes 


L/,(A) = -L/J(A), U:{\*)=U,{\), QW{-\)Q-^ =U,{\). (47) 

For (HSP conditions (I57P are equivalent to 

B = -B'^, B = QB*Q. (48) 


We recall that (see dasp) 

Mi,{u)=In+ 2 + . A- QAQ, A = Qa){a, {a = {i,u^), \u\ = l. (49) 

A — lU \ + IV 


The compatibility condition 

DrMy{u) - Sy{Uy) My{u) + My{u) Uiy = 0 (50) 

for the Lax-Darboux pair 

= M^(ri)4' (51) 

is equivalent to the following two equations 

BiA = AB, (52) 

Ar = Bi-B-^ {BiQAQ + QBiQA - QAQB - AQBQ ), (53) 


where Bi = Si,{B). Indeed, the left hand side of equation (I50p is a rational matrix function of A. 
This function is vanishing at A = oo and conditions ()52p and (|53p are equivalent to the vanishing 
of second and first order poles at A = ±iv. 
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It is easy to check that 

B = —- {Qa-i){a - a){a-iQ ), a_i) = 5”^(a)) 

{aQa^i) 

is a solution of (j52p satisfying the conditions (|48p if 7 is any real function of r. Without any loss 
of generality we can set 7 = 1 (by a point transformation). Substituting the expressions of A and 
B into (1^ . we get 


a i) 


Q 


+ 


(aia) 


{aQa-i) 2{aiQa) 

(oi , (aia) 


Qa) - 


ai) 


(aa i 


+ 


{aiQa) 2(aQo_i 

(a_i (aa_i 


{aiQa) 2{aiQa) 


iaQ + 


{aQa-i) 2{aQa- 


Qa) ) {a 
{aQ 


= 0 . 


Hence, there exists a scalar p such that 


O-t) + 


(a^ - 


a i) 


+ 


(aia) 


(aQa_i) 2{aiQa) 

,Q , 

{aiQa) 2(aiQa) 


. a_i) (aia) . 

{aQa-i) 2{aiQa) ^ \am/ 

Due to (o = the above equation becomes into 


ai) 

(aa i) 

(aiQa) 

2{aQa_i) 

(a i 

(aa i) 

{aQa-i) 

2{aQa_i) 

0 and 


ai) 

, Qa 


Qa) = pa); 
{aQ = -p{a. 


+ 


U-l 




uTu — 1 

+ - -u - 


Ui 


— 1 


u_i + 1 2{u[u + l) u^u + 1 2{u^U-i + l) 


= 0 . 


So the unit vector u satisfies the equation 


_O^C u_i+u 

Uf 2 {Si/ 1) I rp 

\u-^ + U-i 


2 • 


(54) 


We sum up the result in the following theorem: 


Theorem 2. The evolutionary vector equation {54-^ possesses a Lax-Darboux representation i51\) 
with 

Mi/{u) = In+ 2 + ' Qa){a- . a){aQ, {a = {i,u^), u € S"", (55) 

X — iv X + IV 

^ iv Qa-i){a — a){a-iQ iv a-i){aQ — Qa){a-i 
^ X — iv (aQo-i) X + iv (aQa-i) ’ 

where In +2 is the identity matrix of size n + 2, Q = diag{—l, 1,..., 1) and v gM. is a constant. 

Equation (1541) is is known as the discrete Heisenberg spin chain whose link with discrete 
geometry was studied in m- 

It follows from (I35p that the unit vector u is related to the unit vector v by 


Vi+ V 
Ul + u 



In the following proposition, we give the evolutionary equation for the unit vector v. 


(56) 
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(57) 


Proposition 5. Equation (|5^[ ) and equation 

|t 7 _i + v\‘^{vi + 1)) — \vi + v\‘^(v-i + v) 

1 ) ^ - 

+ r>) + |t7_i + + v| 

are related by the Miura transformation 
Proof. Using we can rewrite dSZI) into 

|r)_i + v\^\vi + v\u — It^i + v\^\v-i + v\u-i |d_i + — lui + v\u-i 

'D — - — -^ 

^ (n^u_i + l)|r;_i + t;||Di + r)| vT'u-i + l 

This leads to 

Vir+Vr Ui \V2+Vi\u \v + v\u U_i 

|t7i+rj| vT'ui + l (ri^ui + l)|t;i + t;| (u^u_i + l)|t;i + t;| 1 + 1 


Substituting this into and using the fact that u is the unit vector, we have 


Vir + Vr j,Vir+Vr Ui + U U + U_i U_i + U 

Wt = -:- U - - rU = — -r--r = 2(Si, - l)-r^ - 

I«. _z 


ui + v\ 


Ui + v\ 


Ui + \ U_1 + 1 


ki-' + u_i 


and thus we proved the statement. 


□ 


Both equations (|34l) and (I57p are obtained from the same Darboux matrix My. Thus they share 
the same generalised symmetries and conserved densities derived from the zero curvature conditions 
[la iisi da E], Therefore, equation dSZD can be viewed as a symmetry of the nonevolutionary 
equations (fMl) and (l38l) . 

In 2008, Adler presented the classification of isotropic integrable Volterra-type lattices on the 
sphere m- The author gave a list including six integrable equations denoted by V1-V6. Equation 
(1541) is the V6 when <5 = 1 in the list. Equation (|57p after scaling of r is the V5 when e = 1 and 
k = 2. Here we established the relation between these two equations. 

As pointed out by Adler equation (1541) is integrable without the constraint vector u being 
the unit vector. Introducing the transformation 


U-i + u 
u_l + uP ’ 


we have 


Wr = 2{w'^w){wi — m_i) — Aw'^{wi — (58) 

which is one of vector modified Volterra lattices recently discussed in m- 

Equation ([55D is well studied in m, as an example of integrable Jordan triple. It is a Hamil¬ 
tonian system with the Hamiltonian operator is 

T~Lw = {w^w — 2 wuP"){S — S~^){w'^ w — 2 wnP"). 

and the Hamiltonian function is pw = — In |m|. The authors also gave its master symmetry 

Tw = (n — 1 )wt- + 2{w^ w)wi — wi)w 
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to generate infinitely many commuting generalised symmetries and conserved densities starting 
from Wr and pw 

Using the Miura transformations, we are able to write down the corresponding Hamiltonian 
operators [25], Hamiltonian functions and master symmetries for equations (1541) and (1571) as follows: 


Tiu = {S - 1){S+ 1) \ pu = ln{\ui + u\); 

Ui + u 

Tu — nUj- + 1 ■ pr , 

|rti + up 


Ty = nVr] 


/ (vf + u'^)(u_i + v) 
V \v_l + u||ui + v\ 



where we use the lower-index u and v to indicate their correspondences to equations (I54p and (1571) 
respectively. Using the above master symmetry, we are able to compute higher order symmetries for 
(I57p sharing with both (|34p and (|38p . Recently a new method for construction of master symmetries 
of homogeneous integrable evolution equations (the O-scheme) was proposed in [23] . It would be 
very useful to extend the O-scheme to the classes of equations studied in this paper. 

For the scalar sine-Gordon equation (1421) . the local flow (1571) becomes 


cos ■ 




■ sm 


6>i+6> 


— COS 


6>i-e 


sm 


8 - 1+0 


COS 9 


Oi—O—i I -I 

cos o + 1 


= tan ■ 


-9-1 


which appeared in [23] as a differential-difference version of the modified Korteweg-de Vries equa¬ 
tion. 


5 Yang-Baxter map, integrable partial difference systems and the 
problem of re-factorisation 

It is well known that the problem of re-factorisation of a product of Darboux matrices can be 
associated with the construction of Yang-Baxter maps jSSmsiET!. A Darboux (or discrete Lax) 
representation for integrable partial difference equations can also be seen as the same problem of 
re-factorisation. 

5.1 The Yang-Baxter map associated with the vector sine-Gordon equation. 

Let us consider a product of two Darboux matrices and the problem of re-factorisation 

M,(x)M^(y) = M^(Y)M,(X), p ^ ±iy, x,y,X,YE5’^, (59) 

where My[u) is a Darboux matrix (I55p defined in Theorem [2l On the left hand side of (I59p the first 
and second factors are a matrices with poles at ±i/ and ±/U respectively. On the right hand side the 
order of the factors is opposite. If this problem of a re-factorisation has a unique solution, namely 
for given x,y one can find vectors X, Y uniquely, then it defines a map : S'” x S'"" i->- 5” x 

: (x,y) 1-^ (X(x,y;i/,/r),Y(x,y;i/,^)) , x,y € 5”. (60) 
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Starting from the map we define maps : S'"" x S"" x 5” i-> S'” x S"" x 5”, i < j as follows: 


y^l2 

D,fl 

(x,y,z) 


(X(x, y; v,p),Y (x, y; v, p),z) 

y^13 

(x,y,z) 


(X(x, z; u, K),y,Y (x, z; u, k)), 

y^23 

(x,y,z) 


(x, X(y, z; p,k),Y (y, z; p, k)) 


Definition 1. The map Yy^f^ is called a Yang-Baxter map, if it satisfies the Yang-Baxter relation 

^13 „ \^23 _ v23 „ -v^l3 „ -v^l2 (61) 


o o y23 ^ y23 yl3 yl2 

V,fiL V,K /i,K II,K V,K 


It is easy to show that the map (l60]l associated with the problem of re-factorisation (f5^ satisfies 
the Yang-Baxter relation if the factorisation of the triple product My['x)M^{y)Mf^{z) is unique 


Theorem 3. Suppose that 

(i) the re-factorisation problem f5!9l) defines a unique map Yy^^^; 

(ii) for three different real numbers p,, v, k satisfying p —v, p —k and k —v it follows from 

M,(x)M^(y)M,(z) = M,(X)M^(Y)M4Z) 
that x = X, y = Y, z = Z, that is, the factorisation is unique. 

Then the corresponding map Yy^^ is Yang-Baxter. 

Proof. For any x, y, z € S” it follows from (15911 that 

M,(x)M^(y)M,(z) = M^{Y)My{X)M,{z) = M^(Y)M,(Z)M^(X) = MfiZ)M^{Y)My{X) 


y.t 


(X,Y,Z) 


(X,Y,Z) 


(x,y,z) (X,Y,z) 

Swapping the matrices in a different order we get 

M,(x)M^(y)M«(z) = M4x)M,(Z)M^(Y) = M,(Z)M,(X)M^(Y) = M,(Z)M^(Y)M,(X) 

y23 

(x,y,z) ^ (x,Y,Z) (X,Y,Z) ^ (X,Y,Z) 


Thus we have 

^23 „ -1^13 


■■ (x,y,z)^(X,Y,Z), Yy\^^oYy\ioY^% : (x, y, z) ^ (X, Y, Z). 

The uniqueness of the factorisation that (X, Y,Z) = (X,Y,Z) leads to the required Yang-Baxter 
relation (|UT]) and thus we complete the proof. □ 


Let us compute the Yang-Baxter map (|6UI) corresponding to the Darboux matrix 
Theorem 4. The Yang-Baxter map (Eg) corresponding to the Darboux matrix (Eg) is of the form 

{u^-p^)^ + 2p{u + p{^^y))y 


X{x,y,n,p) = 
Y{x,y;u,p) = 


+ p‘^ + 2pv{xYy) 

{p^ - n^)y + 2v{p + z/(x^y))x 
1/2 + /i2 + 2pv{x^y) 


(62) 

(63) 


where p ±i/ are constant and x,y, X, Y € S”. The map is involutive Yy^^ o Yy^^j^ = id and it has 
the invariant x^y = X^Y. 
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Proof. We derive the map using the re-factorisation of two Darboux matrices (I59p and adapting 
the notations 


a) = 


i 

X 


h) = 


i 

Y 


a) = 


h) = 


At A = oo equation (fM|) is satisfied. Taking the residue of aX \ = iv we get 


Qa){a[l + -^Qh)(b--i^b){hQ ) = ( 1 + 


V — IX 


which implies 


V + IX 


Qa) — ( T 


u — fX 


-Qb){b- 


IX + fx 


^ -Qb){b- -^b){bQ]Qa)T 


{a = T{d{l + 


IX — fX 

fJ- 


IX + fX 


-Qb){b- 




IX — n ' ' ' IX + fx 
for some scalar T 7 ^ 0. Taking the first entry in (1661) we find that 


b){bQ 


T = 


1/2 -^2 


b){bQ]Qa){a, (64) 

(65) 

( 66 ) 

(67) 


1/2 /i2 + 2fxix{-x^y) ■ 

Substitution of (|67)l in (f 66 |) leads to (1621) . Similarly, the residue of (|M|1 at X = ijx yields 
By direct calculation, we can check that x^y = X^Y using the fact x, y € S'”. Indeed, 

+ ^^(x^y)))(x 2 -y) ^ 

(z/2 + ix'^ + 2^i/(x^y))2 ^ 

To prove the involutivity of the map, we simply write the map into the matrix form 


(X, Y) = (x, y) P(x, y), (68) 

where 

p(^ ^ _ I _ ( 2^1^ + 2z/2(x^y) 

’ z/2 -h /i2 + 2^z/(x^y) V 2z/Ai + 2/^2(x'^y) - z/2 

Notice that P(x, y) = P(X,Y) due to the invariance x^y = X^Y. It can be easily checked that 
p2(x, y) = I 2 , which immediately leads to the involutivity of the map. 

It follows from (l68]l that 

(x,y) = (X,Y)P(X,Y), (69) 

which compatible to identities ([651) and the one obtained from the residue of (l5^ at A = fz/ and 
A = ijjL, respectively. Substitution of either ([68]) or (f6^ in (f6^ leads to the identity. 

Finally, we show that the factorisation is unique. We know that the Darboux matrix is invertible 
and (M^(x))-i = Mj(x). Thus if M^(x)M^(y)M«(z) = M^(X)M^(Y)M«(Z), then 

M^(y)M,(z) = Mj (x)M,(X)M^(Y)T4(Z). 

The left hand side is regular at A = z/. The second order pole at A = z/ in the right hand side 
vanishes only if x^X = 1. Since both x and X are unit vectors, we have x = X. Then the right 
hand side of the equation is regular at A = z/. In the same way, we can prove that y = Y and 
z = Z. According to Theorem [3l the map defined by (|62p and (I63p is a Yang-Baxter map. □ 
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Proposition 6. Yang-Baxter map l6^) . [6^) has a linear vector invariant 

vyi + = v'K + nY. (70) 

Proof. The identity can be proved by a direct substitution of (l62]) and ([63]) in (I70]) . Another way 
to prove this statement is to consider the Laurent expansion of equation (|59ll in A at infinity. Then 
(1701) follows immediately from the coefficients at A~^. □ 


5.2 Integrable partial difference system 


Let us consider two Darboux transformations for the vector sine-Gordon system m corresponding 
to Darboux matrices My{^) and M^(Y) with distinct positions of the poles p, / According to 
Theorem [1] the vectors X, y can be expressed in terms of the original variables v = (/3,q:^)^ as 


X = 


Sy{v) +v 


S^{v) + v 


V = ^ (71) 

|5,(n)+n|’ ^ |5^(^) + ^|’ ^ ^ 

where S^, and 5^ are the corresponding shift automorphisms. The Bianchi permutability condition 
for these two Darboux transformations (I16|) has the form 

M,(5^(X))M^(y) = M^(5,(y))M,(X), (72) 

which coincides with (j59|) where 


X = 5^(X) = 


^^ 1,1 +i?o,i 


1 ^ 1,1 +^’ 0,1 

Here we used the notation SfS'ffv = Vnm 


Y = S,{y) = 


^’ 1,1 + 1 ^ 1,0 


(73) 


1 ^ 1,1 + 1 ^ 1,01 

fj_ ^ — ^ri,m introduced in Section [2] and convention do,o = v. Having 
made identifications (1711) , (|73|) and using Theorem S] we can show that the Bianchi permutability 
condition ([72]) is equivalent to a single quadrilateral equation for variables Vij € S'” on the two 
dimensional lattice This equation can be written in a few equivalent forms. 

Theorem 5. Let v E S^, u ^ p & M. and 

{u^-p^)^ + 2p{n + p{^^y))y 


f = 


g = 


1/2 + ^2 2/iz/(x'^y) 

{p'^ - i/2)y + 2i'{p + i/(x^y))x 


z/2 + //2 2pn{-K^y) 


F = 


G = 


where 


X = 


^’ 1,1 + •I’O,! 


y = 


U0,1 + V 


(z/2 - ;,2)x + 2p{u + ^(X^Y))Y 

1/2 + /i 2 + 2pv{X^Y) 

{p^ - j/2)Y + 2v{p + J/(X'^Y))X 

z /2 + ^2 + 2pu{XrY) 

+^’ 1,0 


+ ^’o,ir “ |i’o,i + ^’i 

Then the following equations are equivalent 


X = 


|i’i ,0 +t’| 


Y = 




(74) 


(a) 

Vlfl+V 

li^po + fl 

(b) 

'I’l.l + ^’ 1,0 

- - - - - ~ §? 

(c) 

■^ 0,1 +v 

Ki+^r ’ 

(d) 

1 , 1 = 

(e) 

1 ^ 1,0 = -v + 2 f(f^u), 

(f) 

^’ 1,0 = -^’ 1,1 + 2 g(g^ui,i), 

(g) 

■*^ 0,1 = -u + 2G(G^u), 

(h) 

^’O,! = -^> 1,1 + 2 F(F^ui 4 ) 

hem is 

equivalent to the Bianchi permutability condition ([7^. 
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V + 2f(f^i;). Indeed, it follows 


Proof. Let us show that equation (a) is equivalent to (e): vi^ = — 
from (a) that 


2 P u = 2 




-V = 2- 


vJqV + 1 


= + ^>1 


1 ^ 1 , 0+'*^1 |^’i,o + '^| 

and in particular > 0 and thus (a)=^(e). Since > 0 we have (e)=^>(a): 


^> 1,0 + V _ 2 f(f'^u) _ 
|upo + u| |2f(f'^u)| 


In a similar way one can show that (b)<^=^(f), (c)<;=^(g) and (d)<^=^(h). To demonstrate (e)<^=^(f) 
we show that 

-u + 2f(f'^u) =-upi + 2g(g'^upi) (75) 

by a direct computation. Let |ui^i + uo,i| = Xg and |uo,i + u| = Then 

= + yl = 2{vl^^v+ 1), x^y = + + + (76) 

which implies 

{vl^i+v^)vi,i = Xsys^^y-'^yl, (^^^i + u^i)u = x,y,x^y - (77) 

Using (1761) . (1771) we get 

(zy 2 + + 2yv{y.^y)f{^{f'v) - g(g^i7i,i)) 

= + vl^i)vxf^ + /x(i/ + y.{yFy))yg) {u‘^ - y‘^)xf\vi^i + Uop) 

-2 {{ff^ - u‘^){vq i + v^)vi^iyf^ + i/(/i + z^(x'^y))x5) u{fi + z/(x^y))x7^(ui,i + Uo,i) 

+2 ((z/2 - y‘^)ivJi + u^i)ux7^ + ix{v + y(x^y))y^) y{v + y(x'^y))y7^(uo,i + v) 

- + ^{'y^^y))xs) - i^^)y7H^o,i + v) 

= \{v‘^ + P^ + 2^i/(x^y))2(u - xii^i). 

which leads to the required identity ([75]) . Equation ([M]) is equivalent to (|62]1 . (I63)l . which after 
the identification (fTTjl . (173)) implies that the Bianchi permutability condition ([72]) is equivalent to 
equations (a) and (b). Thus (172)) <^=^(a)<^=^(b)<^=^(e)<;=>(f). In the same way we can show that 
equation ()^ <;=^(c)-^=^(d)-^=^(g)<;=>(h). Observation that (1721) «^=> (l68l) «^=^ (1691) completes the 
proof. □ 


Thus, the condition (1771) is equivalent to a quadrilateral equation (one take any equation from 
the list (a)-(h), depending on the problem) which is a new vector isotropic integrable system with 
dependent variable u G S'” on a sphere. It is natural to refer this system as discrete vector sine- 
Gordon equation. Let us take equation (e) 

1 ^1,0 =+ 2f(f^u) (78) 


as a representative. The Bianchi permutability condition (|72p for two Darboux matrices play the 
role of a Darboux (or discrete Lax) representation for (I78p . 
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Using identification ()7ip . (1731) we can recast the vector invariant (j7np of the Yang-Baxter map 
(Propositional) in a local conservation law for (|78p . that is, 




^1,0 + ^ 


(5. 


+ V 


(79) 


System (j78|) possess an infinite hierarchy of local conservation laws, which can be found using 
formal diagonalisation nanH] of the Darboux matrices or using the master symmetry (15911 . Notice 
that for the equation (1751) the unit vector v satisfies the identity 

(^pi +'»^^i)(^ + '»^o.i) ^ (^pi+^po)(^ + ^i,o) 

1^1,1 + ^>0,ill'll + 1^0,11 + ^’i,o| ’ 


which follows from the invariant x^y = X^Y (Theorem S]). 

For the discrete vector sine-Gordon equation there is a well posed initial value problem with 
initial data given on the staircase 


'^k,k+l I ^ G 'f^k,k T '^k,k+l 7 ^ 0 , lUfc/j| — [u/j— 1 }. 

To find the values vij below the staircase (i > j) one can use equation (e) or (f), for values above 
the staircase (j > i + 1) it is convenient to use equation (g) or (h) and their shifts (by and S^). 
Thus the elimination map [29], which is a useful tool for study symmetries, conservation laws and 
other structures associated with this integrable system can be correctly defined. 

The local symmetry (|57p of the Backlund transformation ()34l) and ()38p results in two symmetries 

+ u) - |ui,o + ^»p(^)l^o + v) 


Dr^V = 


Dr^V = 


(■u^o + vT){v_ifl + u) + |u_i,o + ' 

|^>0 _1 + ^>P(^>o,l + u) - |U04 + ^>p(^>o _1 + v) 
(■u^i + u'^)(uo-1 + u) + |uo-1 + t7||uo,l + v\ 


of the discrete vector sine-Gordon equation. 

Backlund transformations p4p and (I38p give the non-local symmetries of the discrete vector 
sine-Gordon equation as follows: 


Di 


( ^’1.0 + ^’ \ X i X r, f •^0,1+1’ ^ 1 / \ 


T)3,a Qio-Fq: 

[S^ - 1 )—= -2iy- 


Qo.i-Fq: 

- 1)—^ = -2/^7 


/3 " |ui,o /3 |i’o,i 

In the case n = 1, corresponding to the scalar sine-Gordon equation (I42p . we have 


(80) 

(81) 


^>1,0 + V 


Substituting them into (|75|) . we get 


/ 01,0-^6* \ 

/ COS — \ 

/ 

1 ■ 01,0+0 1 ’ 

V 2 / 

x=( 


cos 


sm ■ 




y = 


cos 


sm 


^0,1-1-® 

2 

Op,1+9 


sm 


cos 


01 


o+6> _ sin + 2yu(zv + ^cos 


sm 


^•0,1-1-® 

2 


2 2 I 2 I o 9 i i —9 

COS -+2 — 

/ 2 2\ 9i i-t-00 1 I o / I 9l 1 — 9\ Op 1+0 

0l,p+0 _ ) cos ’ 2 + H COS - ) COS - 


+ IfiVCOS 


OlA-0 
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Using trigonometric identities, we can rewrite it as 


^•1,1+^0,1+^1,0+6 

4 


V cos ' —h ^ cos 


60,1+61,0+30—01,1 


i/sin 


6 i.i+ 0 o,i—01,0 —0 I ■ 00,1+0—01,1—01,o\ _ 


+ /xsin 


= 0 ; 


01 , 1 + 00 , 1 + 01 , 0+6 


V sin ' —h jjL sin 


60,1+61,0+30—01,1 


i/sin 


01,1+00,1-01,0-0 I ■ 00,1+0—01,1—01,o\ 


+ /xsin 


= 0 , 


Since v ^ it follows that 


. 6*1,1 + 6*0,1 — 6*1,0 — 6* . 6*11 + 6*1 0 — 6*0,1 — 6 

V sin-= u sin-, 

4 4 

which is the classical discrete scalar sine-Gordon equation known as the Bianchi-Hirota equation. 


Appendix 


Here we give a proof of the lemma which we use in the proof of Proposition [2j 
Lemma 1. Let a matrix A € Matn+2(C) satisfy 

AA^ = A^A = 0; (82) 

= QAQ, (83) 

where Q = diag{—l, 1,..., 1), then it can be represented as 


H = n^), m, n€ M"’"*'**, |m| = |n| = 1. 

Proof. It follows from ()83p that the form of matrix A is 

e,f€M"+\ H'€Mat„+i(M). (84) 

It follows from ([82]) that 

/i^ = e^e = f'^f; (85) 

HH^ = f f^; H^H = e e^; (86) 

hf-He = he- = 0. (87) 


If /i = 0 then A = 0. Now we assume that /i 7^ 0 and thus e 74 0, f 7^ 0. From (|86p . we obtain that 
rank(i7) = rank(i7FI^) = rank(f f^) = 1 and H can be represented as H = 5fe^, <5 € M. From 
(IHTP we deduce that 6 = h~^. Thus 


A = h 


1 —ih 
ih~^{ h~'^{e^ 



n 


) 


where n = e/h, m = f//i and |n| = |m| = 1 due to ([85p . 


□ 
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